Extended Self- Similarity and Hierarchical Structure in Turbulence 
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It is shown that the two remarkable properties of turbulence, the Extended Self-Similarity (ESS) 
[R. Benzi et ai, Phy. Rev. E 48, R29, (1993)] and the She-Leveque Hierarchical Structure (SLHS) 
[Z.S. She and E. Leveque, Phy. Rev. Lett. 72, 336, (1994)] are related to each other. In particular, 
we have shown that a generalized hierarchical structure together with the most intense structures 
being shock-like give rise to ESS. Our analysis thus suggests that the ESS measured in turbulent flows 
is an indication of the shock-like intense structures. Results of analysis of velocity measurements in 
a pipe-flow turbulence support our conjecture. 

PACS numbers: 47.27.-i 



Fully developed turbulence is characterized by power- 
law dependence of the moments of velocity fluctuations. 
It was suggested by Kolmogorov in 1941 (K41)0 that 
there is a constant rate of energy transfer from large to 
small scales and that the statistical properties of the ve- 
locity difference across a separation r, Sv r = v(x+r)—v(x), 
depend only on the mean energy transfer or equivalently 
the mean energy dissipation rate e and the scale r when 
r is within an inertial range. Dimensional considerations 
then lead to the prediction that the velocity structure 
functions, which are moments of the magnitude of the 
velocity difference, have simple power-law dependence on 
r within the inertial range: 



S p (r) ee (\Sv r \ p ) ~ e p ^ 3 



(1) 



Experiments^) have indicated that there is indeed 
power law scaling in the inertial range but the scaling 
exponents are different fromp/3: 



S p (r) ~ r c " 



(2) 



where £ p has a nonlinear dependence on p. Such a devi- 
ation implies that the functional form of the probability 
density function (pdf) of Sv r depends on r, that is, the 
velocity fluctuations have scale-dependent statistics. Un- 
derstanding this deviation from K41 is essential to our 
fundamental understanding of the small scale statistical 
properties of turbulence. 

Recently, Benzi et aZ.|S] have discovered a remarkable 
new scaling property: S p (r) has a power-law dependence 
on Ss(r) over a range substantially longer than the scal- 
ing range obtained by plotting S p (r) as a function of r. 
This behavior was named Extended Self- Similarity (ESS); 
its discovery has enabled more accurate determination 
of the scaling exponents £ p , particularly at moderately 
high Reynolds numbers assessible experimentally and nu- 
merically, It was later reported that ESS is invalid for 
anisotropic turbulent flows such as atmospheric bound- 
ary layer and channel flow||, || |). This inspires the 



study of a generalized ESS (GESS) , a scaling behavior of 
the normalized structure functions when plotted against 
each other 0], which is still valid in these anisotropic 
flows. The validity of ESS suggests that the different or- 
der structure functions have the same dependence on r 
when r is near the dissipative range || [l(], Very re- 
cently, Yakhot argued that some mean-field approxima- 
tion of the pressure contributions in the Navier-Stokes 
equation would lead to ESSEg). 

A number of phenomenological models have been pro- 
posed to explain the anomalous scaling exponents £ p . 
Among them, a recent successful one is from She and 
Leveque |l3| with a hypothesis of a hierarchical structure 
(HS). When stated for the velocity structure functions, 
the HS hypothesis reads 



S P +i(r) 



= A 
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S P +i(r) 
S P (r) 



[S [o °Hr)] 



1-/3 



(3) 



Here S^(r) = lim^ S p+1 (r)/S p {r) and < f3 < 1 is 
a constant. This hypothesis (and the similar version for 
the local energy dissipation) was supported by experi- 
mental velocity measurements taken in turbulent jets and 
wakes Q |l5| , p~6| . Note that the boundedness of the ve- 
locity will ensure the existence of S(°°'(r). In fact, one 
can show that S^°°\r) = lim p ^ 00 S'p^, and is equal to 
|<5iv| max , the maximum magnitude of (Hv|p"7f. 

It was later reported that the passive scalar structure 
functions |Tg|] and the local passive scalar dissipation Jl9[, 
the temperature structure functions pfj|] and the lo- 
cal temperature dissipation ]2l[] in turbulent convection, 
and the velocity structure functions in a class of shell 
models J22|, Q were all found to possess similar hier- 
archical structures. 

In this Letter, we shall show that the She-Leveque hier- 
archical structure (SLHS) leads to GESS. In other words, 
the SLHS is a special form of GESS. We then give a gen- 
eralized form of SLHS, which is equivalent to GESS, with 
the constant (3 replaced by a function of p. Furthermore, 



2 



it is shown that if |<5iv| max is independent of r for r within 
the GESS range (which will be the case when the most 
intense structures are shock-like), the generalized SLHS 
(and hence GESS) will give rise to ESS. We then conjec- 
ture that the observed ESS in turbulent flows is an indica- 
tion of the most intense structures being shock-like. Con- 
sequently, we predict that in the anisotropic flows where 
GESS but not ESS holds, |<5iv| max has a dependence on 
r. Finally, we present a systematic procedure of analy- 
sis of experimental turbulent signals. Application of this 
analysis to turbulent velocity fluctuations in a pipe flow 
demonstrates that |<5iv| max indeed depends on r in the 
near-wall strong-shear regions where only GESS but not 
ESS is valid but is consistent with being r-independent 
in the centerline of the pipe where ESS is valid. 
First, rewrite Eq. (|h as 



Sk+i(r) 



S fc (r)S(«0(r) 
for integer k, which implies 

S k (r) 
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Si(r) 



S(°°>(r) 



ultiply Eq. (H) for k = 1, 2, . . . , n gives 



S n (r) = B n 



Si(r) 



5(°°)(r) 



[5(°°)(r)] r 



(4) 



(5) 



(6) 



where B n =nj? =1 lri_Q Aj , which further gives 



S n (r) _ B n 



S 3 (r) n / 3 B ™/ 3 



Si(r) 



£M(r) 



3(l-3")-n(l-ff J ) 



Equation (Q) then implies the GESS property: 
T n (r) ~ T m ( r y^ 



(7) 



(8) 



which is a scaling behavior for the normalized structure 
functions, T n (r) = S n (r) / S^fr) 71 / 3 , with the normalized 
exponents p(n, m) even when S n (r) does not have a scal- 
ing behavior in r. For SLHS, 



p(n, m) 



3(1 - (3 n ) - n(l - /3 3 ) 
3(1 - (i m ) -m(l - (3 3 ) 



(9) 



We have thus proved that SLHS implies GESS. 

Earlier, Benzi et al.|| showed that GESS holds if the 
structure functions arc of the form 



S P (r)~ fll (r)^(r) H W 



(10) 



for any functions gi(r), g2(f), and H(p). On the other 
hand, if Eq. (^) with any p(n, m) is valid, we can always 
write the structure functions in the form of Eq. (10) 



say, 



with gi {r) = 5 3 (r)V3, g 2 ( r ) = S q .(r)/S 3 (r)^, 



and H(p) = p(p,q*) for any chosen value of q*. Thus, 
Eq. (£§ is equivalent to GESS. 

Our demonstration here gives the functions g± (r) , g2 (r) 
a meaning. Although the choice of g± (r) and g2 (r) is not 
unique, the following form can be obtained from Eq. (^): 
9l (r) = S , (°°)(r) and g 2 (r) = S 3 {r)/[S^(r)] 3 . Here, 
gi(r) describes the r-dependence of very strong fluctu- 
ations S(°°>(r) and 32 ( r ) describes the normalized in- 
dependence of (typical) weak fluctuations [e.g. 5 3 (r)] 
by [S , (°°)(r)] 3 . Express g± and g 2 this way, we have 



S p (r)~ [S<~>(r) 



5s(r) 



[S(oo)( r )]J 



f(p) 



(11) 



where /(p) is a function to be discussed below. Note that 
when c^M is constant, or the weak fluctuations have the 
same r-dependence as the strong fluctuations, we have 
the K41 scaling. The function f(p) has to satisfy vari- 
ous conditions. By definition, /(0) = 0, /(3) = 1, and 
limp^oo f(p + 1) - f(p) = linip^oo f(p)/p = 0. Fur- 
thermore, the boundedness of the velocity restricts that 
df(p)/dp > 0. 

Note that Eq. (|Tl|) is a general expression of GESS. 
Rewriting Eq. (O) in the form of Eq. (0) , we have 



Sp+i(r) 



V)~ Cp 



S P (r) 



S p -i{r)S(°°)(r) 



a(p) 



(12) 



where g(p) = [f(p+l)-f(p)]/[f(p)-f(p-l)]. It is clear 
that SLHS corresponds to the particular case of g(p) = (3 
or equivalently 



f(p) = 



1 - (3 p 
1-P 3 



(13) 



We can list several additional possible classes of f(p), 
which include f(p) = [(pa + l) a - l]/[(3cr + 1)° - 1] and 
f(p) = \n(pa + l)/ln(3cr + 1) with a > and < a < 
1. The latter is the limit of the former when a — > 0. 
These two cases correspond to the two results derived by 
Novikov, Eqs. (18) and (19) in Ref . p(| , using the theory 
of infinitely divisible distributions. 

To study when GESS would further give rise to ESS, 
we rewrite Eq. (pi) in the following form: 



S P (r) 



5 3 (r)V3 



S 3 (r)V3 



with 



Pip, ?) = 



/(p)-p/3 
/(«) " 9/3 



(14) 



(15) 



Thus we see that if there exists an p* / such that 
Sp*(r) 1 ' p is independent of r, then we would have a 
scaling behavior of S p vs S3 (and thus S q for q p). If 
p* is finite, £ p » = which implies Cm = f° r < m < p* 
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because has to be a non-decreasing function of n in or- 
der for the velocity field to be bounded Q . It further gives 
C p = for all values of p if Q p is an analytic function of 
p. This could be avoided only if p* approaches oo, which 
means that S {cc \r) is independent of r. Hence, GESS 
together with the condition that 5(°°)(r) or equivalently 
\Sv r \ maK being independent of r would give rise to ESS: 



S P (r) ~ S 3 (r)"<* 3 > 



(16) 



with ri(p, 3) = f(p) even when S p (r) does not have a 
power-law dependence on r. We note that the pres- 
ence of shocks with maximum velocity, i> ma x, in oppo- 
site directions across the shock discontinuity would give 
|(Ji; r | max = 2ii max , which is thus independent of r. 
Using Eq. (O), an independence of S^°°^ on r implies 



C P = Cs/(p) 



(17) 



For SLHS, linip^ /(p) = l/(l - /3 3 ) hence Eq. (|l7|) fur- 
ther implies a saturation of £ p as p — > oc. A connection 
of the saturation of the exponents with the existence of 
shocks was suggested earlier by Chen and Caop7|. 

Both ESS and SLHS have been observed in a variety 
of fluctuating flow fields. The above analysis points to a 
possibility that the two combined may be indicative of a 
property of the flow field: S^°°'(r) is independent of r. 
One way to check the plausibility of these ideas is to ex- 
amine experimental velocity measurements that demon- 
strate GESS, and to investigate whether |<5« T .| max is in- 
deed independent of r when ESS is valid and otherwise 
dependent on r when ESS is not valid. Note, however, 
that the detectable |(5u r | max in any finite measurement 
would almost surely underestimate the true value. 

If the GESS is of the special form of SLHS, a sys- 
tematic procedure can be developed to perform an indi- 
rect estimate of the r-dependence of 5(°°) (r) as follows. 
First, one verifies if the SLHS is valid by performing a 
/3-test. It consists of computing the normalized structure 
functions T p (r) and obtaining the normalized exponents 
p(j>,q) by measuring the slopes of ln(T p ) vs ln(T g ). Let 
Ap(p, q) = p{p + 1, q) — p{p, q). It is easy to derive the 
following equation when SLHS is valid: 



Ap(p + l,q)=(3Ap(p,q)- 



q(l-0*)-3(l-0*) 



(18) 



If one finds parallel straight lines when ploting Ap{p + 
l,q) vs Ap(p,q) for a set of values of q, we say that 
the data passes the /?— test and the turbulent flow field 
possesses the SLHS property. The slope and intercept 
provides a double estimate of the constant (3. With the 



estimated /3, one can then construct f(p) using Eq. (13) 
For an indirect estimate of S^°°\r), we introduce 



F p (r ru) - HS P (r)/S p (r )]-f(j>)hi[S 3 (r)/S 3 (r )] 



p-3/(p) 



From Eq. fllip, F p {r, ro) should be independent of p and 
equal to log[SM( r )/S(«0( ro )] for r within the GESS 
range. Thus, one computes F p (r, r ) for a fixed value of 
ro, within the GESS range, and plots it as a function 
of r for a set of values of p to estimate indirectly the 
r-dependence of InlS^ (r) / (r )]- In particular, if 
g(oo) - 1S independent of r, then F p (r,ro) = for r within 
the GESS range. 

We have applied the above procedure to analyze hot- 
wire measurements of longitudinal velocity fluctuations 
in a pipe flow |p8f . The pipe is 22.5 m long with an in- 
ner diameter of 10.5 cm, and the Reynolds number is 
1.35xl0 5 . The velocity measurements studied were taken 
at a cross section at 18.2 m away from the entrance both 
at the centerline of the pipe and at a distance of 0.1 mm 
from the pipe-wall. The number of measurements taken 
at each point is 5.76 x 10 7 and the maximum order of 
moments computed is p = 10. We have checked that 
GESS is valid for both sets of measurements but ESS is 
valid only at the centerline. The range of GESS (and also 
ESS when valid) is r « 10 — 500, in units of the sampling 
time 1/48 ms. All r's quoted below will be in the same 
units. We have obtained p(p, q) for the two locations 
and r](j), 3) for the centerline location. In Fig. [l], we plot 
Ap(p +l,q) vs Ap(p, q) for some values of q for both lo- 
cations. The data can be fitted by parallel straight lines 
showing that the hierarchical structure is indeed of the 
SL form in both locations. We estimate the value of (3 
simultaneously from the slope and the intercept and get 
0.93 ± 0.01 and 0.85 ± 0.01 respectively for the center- 
line and the near- wall location. In the inset, we see good 
agreement of r](p, 3) with f(p) for the centerline mea- 
surements. This implies the r-independence of S^-°°\ in 
accord with the validity of ESS in these measurements. 

We have next computed F p (r, ro) with ro = 69 for both 
locations and found that the data indeed collapse when 
p is large. In Fig. @, we plot F p (r,r ) as a function of 
r. We see that F p (r,ro) for the centerline measurements 
is almost zero for r > 30 while that for the near-wall 
measurements shows a clear r-dependence throughout. 
Since ESS is valid for the centerline but not the near- 
wall measurements (see the inset), these results support 
our conjecture that the validity of ESS in turbulent flows 
is an indication of S^ 00 ' being independent of r. 

In summary, our analysis proposes a link between a 
measurable statistical property of turbulence, the ESS, 
and a property of the most intense structure, the r- 
independence of |i5w r | lnax or S^^r). We have devel- 
oped a systematic procedure to test this conjecture when 
the SL hierarchical structure holds. The analyis of ex- 
perimental pipe flow data support our conjecture: the 
near-wall strong-shear turbulence contains more complex 
structures while the centerline fully developed turbulence 
has shock-like structure for S^°°'(r). A consequence of 
the r-independence of S^°°'(r) is the saturation of the 
scaling exponents of the velocity structure functions at 
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FIG. 1: Ap(p + l,q) vs Ap(p,q) for q = l (squares), q = 2 (cir- 
cles), q — 4 (triangles), q = 5 (plusses), and q = S (diamonds). 
Above the solid line (slope = 0.93) are the centerline data and 
below the dashed line (slope = 0.85) are the near-wall data 
(shifted for clarity). In the inset, we compare r/(p, 3) (circles) 
with f(p) (solid line) for the centerline data. 
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FIG. 2: F p (r, ro) for p > 4 for both the centerline and near- 
wall measurements. In the inset, we show the local slope of 
In S^r) vs In £3(7-) as a function of In 5*3 (r), illustrating that 
ESS is valid for the centerline (solid line) but not the near-wall 
(dashed line) measurements. 

very high orders. Further experimental and numerical 
tests are highly desirable. 
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